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\^ Abstract 

,—1 We investigate the effect of supersymmetric (SUSY) threshold corrections on the values 

1-^ of the running quark and charged lepton masses at the GUT scale within the large tan /3 

SpH regime of the MSSM. In addition to the typically dominant SUSY QCD contributions 

pi I for the quarks, we also include the electroweak contributions for quarks and leptons and 

O show that they can have significant effects. We provide the GUT scale ranges of quark 

and charged lepton Yukawa couplings as well as of the ratios to^/tos, nie/rnd, Vr/yb 
and yt/yt for three example ranges of SUSY parameters. We discuss how the enlarged 
04 ranges due to threshold effects might open up new possibilities for constructing GUT 

^^ models of fermion masses and mixings. 
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1 Introduction 

The unification of tlie fundamental forces of the Standard Model (SM) is one of the guiding 
principles in the search for a more fundamental theory of nature. In addition to providing 
a unified origin of the gauge interactions, Grand Unified Theories (GUTs), based e.g. on 
the gauge symmetry groups SU(5) [1] or SO(IO) [2], also unify quarks and leptons of the 
SM in representations of the unified gauge groups. This property makes them attractive 
frameworks to address the flavour problem, i.e. the question of the origin of the observed 
pattern of fermion masses and mixings. Another attractive feature of left-right symmetric 
GUTs is the appearance of right-handed neutrinos in their particle spectra, which become 
massive after spontaneous symmetry breaking to the SM and thereby lead to the small 
observed neutrino masses via the seesaw mechanism [3] . In order to solve the gauge hierarchy 
problem inherent in high-energy extensions of the SM and to make the running gauge 
couplings meet (at the so-called GUT scale Mqut ~ 2 x 10^^ GeV), the idea of Grand 
Unification is typically combined with that of low-energy supersymmetry (SUSY). 

In order to construct successful GUT models of flavour, it is desirable to know the 
approximate GUT scale values of the quark and lepton masses and mixing parameters. 
The experimental data on the masses of the strange quark and the muon, for example, 
extrapolated to the GUT scale by means of renormalisation group (RG) running within 
the SM, gave rise to the so-called Georgi Jarlskog (GJ) relations [4] of m^/ms = 3 and 
ine/iTiii = 1/3 at the GUT scale, which can be realised from a Clebsch-Gordan factor 
after GUT symmetry breaking. The GJ relations have become a popular building block 
in many classes of unified flavour models. In SUSY GUTs, another intriguing possibility 
emerges, which is the unification of all third family Yukawa couplings, i.e. of yt, yt, Vt 
and furthermore, in the context of the seesaw mechanism, of y^y at the GUT scale. To 
investigate whether this relation can be realised in a given model of low-energy SUSY, it 
is well known that a careful inclusion of SUSY threshold corrections is required [5, 6, 7, 8]. 
These threshold effects are particularly relevant in the case of large tan/3, where yt = yb = Vt 
seems achievable. Despite the possible importance of the SUSY threshold effects, in studies 
which interpolate the running fermion masses to the GUT scale these effects are typically 
ignored (see, e.g., [9, 10]). 

The effects of SUSY threshold corrections on the possibility of third family Yukawa uni- 
fication yt = yb = Vt (and also on the less restrictive possibility yb = yr which often emerges 
in SU(5) GUTs) has been extensively studied in the literature (see, e.g., [6, 7, 11, 12, 13]). 
Furthermore, recent studies [14] have addressed the phenomenological viability of this rela- 
tion and have pointed out that under certain assumptions on the soft breaking parameters 
at the GUT scale, yt = yt = Ut may be already quite challenged by the experimental data 
from B physics. SUSY threshold effects on the GJ relations have been discussed recently 
in [15]. Taking into account the typically leading contribution from SUSY QCD loops, it 
has been shown that the GJ relation can be realised under certain conditions on the SUSY 
parameters which govern this correction. 

The main purpose of this study is to include SUSY threshold corrections and calculate 
the possible GUT scale ranges for quark and charged lepton Yukawa couplings as well as for 
the ratios m^/ms, me/md, yr/Ub and yt/yb, which are important input parameters for GUT 



model building. Regarding m^/ms, compared to [15] we include additional corrections from 
electroweak loops with binos and winos for quarks and charged leptons, which, as we will 
show, can have significant impact. Furthermore, instead of trying to fit the GJ relations 
by a sparticle spectrum, our aim is to analyse which alternative GUT scale relations may 
be possible and whether the GJ relation lies within the projected GUT scale ranges. We 
also investigate other relevant quantities of interest for GUT model building, such as, for 
example, the ratio of electron mass over down quark mass, me/rrid-, and of the third family 
Yukawa couplings j/t, Vb and yr- 

The remainder of the paper is organised as follows: In section 2 we discuss the various 
tan /3-enhanced contributions to the SUSY threshold corrections and how we implement 
them within the renormalisation group running procedure of the Yukawa couplings. Section 
3 contains our results for the GUT scale ranges for the quark and charged lepton Yukawa 
couplings as well as for the ratios m^/iris, mg/mii, yr/Vb and yt/yb, and in section 4 we 
discuss possible implications for GUT model building. Section 5 concludes the paper. 



2 SUSY Threshold Corrections 

For large tan/? in the MSSM it is well known that certain supersymmetric one-loop cor- 
rections are enhanced (see, e.g., [5, 6, 7]). The Feynman diagrams corresponding to these 
one-loop corrections are shown in figures 1 and 2. Qi denotes the quark doublets, di and m 
the quark singlets, Li the lepton doublets and Cj the lepton singlets of the i-th generation. 
Hu is the up- type Higgs doublet, Hd the down- type Higgs doublet and superpartners are 
marked with a tilde. For instance, G labels the gluino with mass M3, B the bino with mass 
Ml and W the wino with mass M2. 

After matching the SM with the MSSM at the SUSY scale at one-loop [16, 17], we 
obtain for the affected Yukawa couplings of down-type quarks and charged leptons 

MSSM ^ Vi (2 1 "l 

^^ cos/3(l + eitan/3)' ^ ' ^ 

where yf^ = fhi/v is the SM Yukawa coupling related to the MS-mass rfij of the respective 
particle. With Vu and Vd being the vacuum expectation values (vevs) of the up-type and 
down-type Higgs doublets, tan/3 is defined as tan/3 = Vu/vd and the vev v of the SM Higgs 
is given by f^ = v'^ + v'^. The quantities e^ which govern the tan /3-enhanced corrections, 
with i = d, s, b, e, fi and r, will be specified in section 2.1. Note that since the calculation is 
one-loop, we can use MS-quantities as well as DR-quantities. For large tan/3 the correction 
ej tan /3 can become quite large, so higher order calculations seem to be necessary, however 
it can be shown that the tan /3-enhanced vertex corrections are absent in higher orders [16]. 
In the following, we will restrict ourselves to the case of large tan /3. 

For the up-type quarks (and neutrinos) there are no tan /3-enhanced corrections, and 
the remaining corrections are negligibly small. We will therefore match the up-type quark 
Yukawa couplings using the tree-level relation y¥SSM _ y9™i/sin/3, with i = u, c and t. 
For the neutrino sector we will assume that the small masses are generated by the (type I) 
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Figure 1: Feynman diagrams contributing to the SUSY threshold corrections to down-type quark Yukawa 
couphngs. 

seesaw mechanism [3] operating at high energies close to the GUT scale. In the following, 
we will drop the MSSM-label for the Yukawa couplings and imply yi = y¥SSM_ 



2.1 Corrections to Quark and Lepton Yukawa Couplings 

We give now explicit formulae for the quantities e^ in equation (2.1), which in addition to 
tan/3 govern the size of the threshold corrections. Details of the calculation for the quarks 
can be found in [16, 17]. The formulae are valid for the case of unbroken SU{2)l x f7(l)y 
symmetry, which is appropriate for the situation that the sparticle masses are in the TeV 
range and that we use matching conditions at these high energies above the electroweak 
(EW) scale Mew- 

Turning to the corrections for the down-type quarks we can decompose ej = ep + e^ + 
ef + ey6,b, with [18] 
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Figure 2: Feynman diagrams contributing to the SUSY threshold corrections to charged lepton Yukawa 
couplings. 
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myM^, Vj = my^j.^, wj 
i = d,s,b and where ah mass parameters are assumed to be reah The correction e^ is only 
relevant for the 6-quarks because of the strong hierarchy of the quark Yukawa couplings. 
The function H2 is defined as 
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Note that H2 is negative for positive x and y and \H2\ is maximal, if its arguments are 
minimal, and vice versa. 

The corrections for the charged leptons stem from diagrams similar to the ones for the 
quarks and are shown in figure 2. One difference between the corrections for quarks and 
charged leptons is of course that the SUSY QCD loop contributions ef are absent. Another 
difference concerns the contributions ef with binos in the loops, where due to the different 
hypercharge of the charged leptons the prefactors for these contributions are changed. In 
the last term in equation (2.7) this causes an enhancement by a factor of —9 compared to 
the corresponding term in the quark sector in equation (2.3). The contribution from the 



diagrams with winos ej^, on the other hand, is equal for quarks and leptons. A further 
difference between the corrections for quarks and charged leptons is that in the considered 
seesaw framework the r-lepton Yukawa coupling does not have a relevant correction of the 
e^ type because the corresponding vertex correction is suppressed by the heavy mass scale 
of the right-handed neutrinos. For the corrections for the charged leptons we find 
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for i = e, fi, T. 



2.2 Renormalisation Group Running from EW to the GUT Scale 

The calculation of the Yukawa couplings of the quarks and charged leptons at the GUT 
scale can be accomplished by solving the corresponding renormalisation group equations 
(RGEs) in the (type I) seesaw scenario from low-energy to the GUT scale. For this we use 
the package REAP introduced in [19] , where also a summary of the relevant RGEs can be 
found. 

For our analysis, we take as input values the running quark and lepton masses at the 
top scale fht{fht), which have been calculated recently in [10] with up-to-date experimental 
values for the low-energy quark and charged lepton masses, and evolve them first to the 
SUSY scale MsusY using the SM RGEs. At the SUSY scale we match the SM with the 
MSSM to obtain the running MS Yukawa couplings via equation (2.1). Since we consider 
one-loop running, we can neglect issues of scheme dependence such as transformations from 
MS to DR quantities. Two-loop running (and scheme-dependent) effects are small compared 
to the tan /^-enhanced threshold corrections and can be neglected. 

As the next step, we solve the RGEs from the SUSY scale to Mqut taking into ac- 
count possible intermediate right-handed neutrino thresholds as discussed in [20]. For our 
numerical calculations we use REAP, which solves the complete set of one-loop RGEs and 
automatically includes the right-handed neutrino thresholds. We will comment on the pos- 
sible effects of right-handed neutrino thresholds, which depend on the additional degrees of 
freedom in seesaw models, in section 3.4. If not stated otherwise, they are ignored in our 
analysis. 

We note that there are SUSY scenarios which may lead to corrections to our approach 
of one-step matching at the SUSY scale in the EW unbroken phase. For example, if the 
sparticle spectrum is light, effects of EW symmetry breaking may have to be taken into 
account for the calculation of the SUSY threshold corrections. Another example is the 
possibility of having a (relatively) split sparticle spectrum, in which case matching at one 
scale would be a bad approximation. When we present explicit examples in the following, 
we choose parameters where our assumptions are justified to a good approximation. 
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SUSY parameter 


Case (7+ 


Case g- 


Case a 


ruf in TeV 


[0.5, 1.5] 


[0.5, 1.5] 


[0.5, 1.5] 


Ml in TeV 


[0.5, 1] 


[0.5, 1] 


[1.65, 3.3] 


M2 in TeV 


[1,2] 


[1,2] 


[0.5, 1] 


M3 in TeV 


[3,6] 


[3,6] 


[-9, -4.5] 


/i in TeV 


0.5 


-0.5 


0.5 


At in TeV 


±1 


±1 


±1 


^SUSY in TeV 


1 


1 


1 



Table 1: Example ranges of SUSY parameters (at the matching scale Msusy) used in our analyses in 
sections 3.1 and 3.5. The choices of gaugino masses in the cases g± are inspired by universal gaugino masses 
at the GUT scale and in case a by anomaly-mediated SUSY breaking. We therefore use the low-energy 
approximation Mi : M2 : M3 = 1:2:6 for the cases g± and Mi : M2 : M3 = 3.3 : 1 : —9 for case a as a 
constraint. 

3 Quark and Lepton Yukawa Couplings at the GUT Scale 

In this part we present our analytical and numerical results. We start with a semianalytic 
discussion of the individual contributions to the threshold corrections presented in equations 
(2.2) to (2.8), which allows to estimate their size as well as the size of the total corrections 
ej. We will then turn to the numerical analysis, where we quantitatively discuss the effect 
of the most relevant parameters and present our final results for the quark and charged 
lepton Yukawa couplings (and mass ratios-*^) at the GUT scale. To isolate the effects of the 
SUSY (MSSM) parameters from the uncertainties induced by the quark mass errors, we 
first use best-fit values for the low-energy fermion masses. These errors are later included 
in our final results presented in tables 6 and 7. We note that the quark mass errors have a 
significant effect, whereas the charged lepton mass errors are negligibly small. 



3.1 Semianalytical Approach 

For our semianalytical approach we first define three example choices of possible ranges of 
the relevant SUSY breaking parameters. They are listed in table 1 and will be referred 
to as cases g+, g- and a. The cases (7+ and g- are inspired by scenarios with universal 
boundary conditions for the gauginos (where Mi/gf = M2/g2 = M^/g'^) and case a by 
anomaly-mediated SUSY breaking [21, 22] (where Mi/{gf/3i) = M2/{gl(32) = M^/{gll3^), 
with (/?!, /32,/33) = (33/5,1,-3)). We note that instead of these relations at the SUSY 
scale, we use the low-energy approximation Mi : M2 : M3 ~ 1 : 2 : 6 for the cases g± and 
Ml : M2 : M3 ~ 3.3 : 1 : —9 for case a. We furthermore only introduce relations for the 
gaugino masses, whereas for the sfermion masses m? and the /i and At parameters we do not 
apply restrictions from a specific model of SUSY breaking. The example ranges are chosen 
such that we can neglect effects which are suppressed by Mew/^^SUSY, such as mixing 



^We note that when we refer to fermion masses at the GUT scale, what we mean is simply the Yukawa 
coupling multiplied by the low-energy value of the corresponding Higgs vev. 





Case 5f+ 


Case g^ 


Case a 


e^ in 10"-^ 


[3.52, 9.31] 


[-9.31, -3.52] 


[-7.85, -3.10] 


ef for quarks in 10^"^ 


[-0.31, -0.08] 


[0.08, 0.31] 


[-0.30, -0.16] 


ef for leptons in 10"^ 


[-0.18, 0.30] 


[-0.30, 0.18] 


[0.01, 0.30] 


ef in 10-^ 


[-2.21, -0.98] 


[0.98, 2.21] 


[-2.21, -0.72] 


ey in 10-3 


sign(^t) [0.84, 4.65] 


-sign(Ai) [0.84, 4.65] 


sign(^t) [0.84, 4.65] 



Table 2: Ranges for the various contributions to the SUSY threshold corrections (corresponding to the 
example ranges of SUSY parameters in table 1) for i = d,s,b and i = e,fi, r, respectively. For the charged 
leptons there is no contribution ef and e". The ranges for ef' are the same for quarks and leptons. 



between left- and right-handed sfermions, and that our approach of one-step matching is 
justified to a good approximation. We note that left-right mixing effects may nevertheless be 
important [17], for instance in so-called inverted scalar mass hierarchy (ISMH) scenarios. We 
also note that without specifying the remaining SUSY parameters (which do not enter the 
formulae for the threshold corrections) we cannot apply various relevant phenomenological 
constraints on the spectrum. To do so, we would have to impose further constraints on the 
soft breaking parameters (as done e.g. in [14]) which is, however, not our intention. 

For the parameter ranges specified in table 1, we can estimate the corresponding ranges 
for the different components of e^. Scanning over the ranges of SUSY parameters in table 1 
we obtain the resulting ranges presented in table 2. The gauge couplings are evaluated at 
the SUSY scale. From these ranges we can already see that the electroweak contributions 
cannot be neglected. In our scan we find that they can amount up to about 50 % of the 
QCD contribution, larger than one might suspect from table 1 due to correlations between 
the corrections. We can also see from table 2 that, for the quarks, in case a the inclusion 
of the electroweak corrections results in an enlargement of the threshold correction because 
here the QCD and the electroweak corrections add up, whereas in case g± it results in a 
reduction of the total correction because both contributions partially cancel. 

The ranges for the Cj from table 2 can now be used to obtain (naive) analytic estimates 
for the ratios of the fermion masses (or Yukawa couplings) at the GUT scale. For example, 
in leading order the GUT scale ratio me/rrid is given by 

me(MGUT) me(MGUT) 1 + edtan/3 me(MGUT) 



mrf(MGUT) md(MGUT) 1 + eetan/3 mrf(MGUT) 



(1 + (ed - €e) tan (3) + 0{el tan^ /3), 

(3.9) 

where m(MGUT) denotes the fermion masses at the GUT scale without SUSY thresholds. 
We will use the analogous formula for the second generation. For the third generation we 
take the ratios of the Yukawa couplings so we have to take into account an additional tan (3 
factor for yt/Ub- We will later on compare these estimates with the numerical results for 
the same ranges of MSSM parameters (cf. table (5)). For the values of the masses and 
Yukawa couplings at the GUT scale, we take the values calculated with REAP setting all 
SUSY threshold corrections to zero and using the best-fit values for the fermion masses as 
low-energy input. These values for the Yukawa couplings are collected in table 3. In the 
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Ve in 10"* 


y^ in lO-"-^ 


Vr 


tan/3 = 30 


0.62 


1.31 


0.23 


tan /? = 40 


0.88 


1.85 


0.34 


tan /? = 50 


1.21 


2.55 


0.51 






Vd in 10-4 


Vs in 10-2 


Vb 


tan/3 = 30 


1.57 


0.30 


0.18 


tan /? = 40 


2.22 


0.43 


0.26 


tan /3 = 50 


3.06 


0.59 


0.39 






Vu in 10-6 


Vc in 10-4 


Vt 


tan /? = 30 


2.73 


1.33 


0.49 


tan /3 = 40 


2.75 


1.34 


0.50 


tan (3 = 50 


2.77 


1.35 


0.52 



Table 3: Best-fit values for the Yukawa couplings at the GUT scale without SUSY threshold corrections 
(case 0) for Msusy = 1 TeV and different values of tan/3. 





Case 


Case (7+ 


Case g^ 


Case a 


rrie/md 


0.39 


[0.35, 0.64] 


[0.15, 0.44] 


[0.16, 0.44] 


m^/nis 


4.35 


[3.83, 7.01] 


[1.69, 4.87] 


[1.81, 4.85] 


Vr/yb 


1.32 


[1.16, 2.13] 


< 1.48 


< 1.47 


yt/vb 


1.93 


[1.65, 2.92] 


< 2.21 


< 1.98 



Table 4: Semi-analytic (naive) estimates for the ranges of the mass and Yukawa coupling ratios at the GUT 
scale (corresponding to the example ranges of SUSY parameters in table 1) for tan/3 = 40. Case refers to 
the case without SUSY threshold corrections. For the ranges involving yi,, the lower boundaries depend on 
the cut we had to introduce in order to keep yt perturbative up to Mgut and therefore have been omitted. 



following (e.g. in table 4) we will refer to the case without SUSY thresholds as case 0. 

The results of these estimates are collected in table 4. We note that these estimates are 
naive in the sense that we have combined the maximal and minimal values of each of the 
contributions to e,, neglecting possible correlations between them. For example, we do not 
account for the effect that the QCD corrections become large, if M3 and thus the gaugino 
masses are large and the sfermion masses small, whereas the wino corrections become large, 
if the gaugino masses and the sfermion masses are small. However, as we will later see, the 
estimates nevertheless work surprisingly well. Another effect which we can immediately see 
from the analytic estimates is that yb can become nonperturbatively large if €b tan /3 is large 
and negative. In fact, it can even occur that ejtan/? < —1, which spoils the perturbative 
expansion. Wherever non-perturbative values of 1/5 occur, we only give the upper boundaries 
of the ranges yr/yb and yt/yt and the lower boundary for yi, itself. 

The naive estimates already suggest that with SUSY thresholds included, a wide range 



of GUT scale values of down-type quark and charged lepton masses (or Yukawa couplings) 
could be realised. Taking a preliminary look at the predicted ratio for m^/ms, the naive 
estimates suggest that with the SUSY parameters of example g+ , the GUT scale value of 
m^/nis is typically significantly larger than the GJ relation of m^/ms = 3. On the other 
hand, scenario g^ and a are well compatible with the GJ relation. Beyond the GJ relation, 
the naive estimates also imply that with SUSY thresholds included, a large variety of GUT 
model predictions for these ratios might be compatible with the low-energy data on quark 
and lepton masses. A full numerical analysis for the example SUSY parameter ranges g± 
and a will be presented in section 3.5. 



3.2 Dependence on ^ and tan/3 

Before we proceed with the numerical analysis for the example ranges of MSSM para- 
meters of table 1, we discuss the dependence on ^ and tan/3, which has been kept fixed 
in the last section. The dependence on ^ is rather important, because all corrections are 
proportional to ^ or 1/^. The parameter ^ therefore gives the overall sign of the corrections 
and determines if the Yukawa couplings are enhanced or reduced by the SUSY threshold 
effects. In addition, tan/3 is very important because the threshold corrections are almost 
linear in tan/3 and also because for successful third family Yukawa unification we need a 
large value of tan/3. To isolate the effects of these parameters, we have effectively turned 
off the right-handed neutrino threshold effects, put At to zero and all the other soft SUSY 
breaking parameters and the SUSY scale to 1 TeV (with both signs allowed for M3 but 
with Afi > 0, M2 > 0). In figures 3 and 4 the numerical results are presented as contour 
plots in the /U-tan/3 plane for the four ratios ra^jrad^ m^/ms^ VrlVh and ytlvh for different 
combinations of the sign of /U and M3 . 

There are several interesting points we would like to remark on: First, the overall 
dependence in the plots illustrates the anticipated behaviour from the fact that the leading 
contribution from gluino loops is proportional to [i and that the overall size of the corrections 
is proportional to tan/3. They also illustrate that for [iM^ < (second and third rows in 
the figures) the total corrections enhance the down-type quark Yukawa couplings leading 
to more stringent restrictions for the possible values of tan /? from perturbativity of yb up 
to the GUT scale. On the other hand, for fiM^ > and M3 > (first row in the figures) 
the total corrections lower the down-type quark masses and in principle larger values for 
tan/3 are possible. Second, interesting conclusions can also be drawn from comparing the 
second to the third row of the figures. From the leading SUSY QCD contribution which is 
invariant under a simultaneous change of sign in fi and M3, one might expect that the plots 
in the second and third rows look very similar (if understood as results for \lJ-\). Differences 
are entirely induced by the contributions from wino and bino loops, since we have chosen 
At = 0. Inspecting the numerical results shows significant differences in the results for ^u < 
and M3 > and ^ > and M3 < 0, which confirms that the EW contributions are indeed 
important and cannot be ignored (as we have already concluded from our semianalytic 
estimates in section 3.1). 
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Figure 3: Contour plots for the GUT scale ratios me/md (left side) and in^/ms (right side) for M3 > 0, 
pi > (first row), A/3 > 0, /i < (second row) and M3 < 0, /.i > (third row) in the ^-tan/3 plane. In the 
grey areas labeled with NP the value of y;, becomes nonperturbatively large. 
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Figure 4: Contour plots for the GUT scale ratios j/r/j/t. (left side) and ytjyb (right side) for Mz > 0, /i > 
(first row), A/3 > 0, ^ < (second row) and M3 < 0, /x > (third row) in the /.i-tan/3 plane. In the grey 
areas labeled with NP the value of j/(, becomes nonperturbatively large. 
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3.3 Dependence on AI^susy and At 

The GUT scale values of the quark and lepton Yukawa couplings also depend on MsusY 
and At- While the correction to the bottom Yukawa coupling can be significant (as can be 
seen from table 2), the effects on the down and strange quark Yukawa couplings are quite 
weak since they only stem from indirect effects (modified RG evolution) due to the change 
of the bottom Yukawa coupling. 

We also looked explicitly at the dependence on MsusY by fixing all other parameters and 
varying only MsusY- We found that changing MsusY can have some effect on the GUT 
scale value of the Yukawa couplings due to the difference in the RGEs between SM and 
MSSM; however this effect is typically much smaller than the uncertainty induced by the 
quark mass errors and the sparticle spectrum, and it nearly cancels out when we consider 
ratios of masses or Yukawa couplings. We have therefore fixed MsusY to 1 TeV in our 
numerical examples. 

3.4 Right-handed Neutrino Threshold Effects 

For analysing the possible dependence on threshold effects from the right-handed neutrino 
sector, we have taken the three examples for SUSY parameters from our analysis on the fi 
and tan/3 dependence, but fixed fi to ±0.5 TeV and tan/3 to 40. For these example param- 
eter points we have investigated the effects for the three different scenarios of sequential 
dominance [23] also used as examples in [24] . With largest neutrino Yukawa couplings being 
0(1), we found that the deviations are typically smaller than 5 %, which is small compared 
to the SUSY threshold effects (in our examples) and also compared to the uncertainties 
induced by the present quark mass errors, especially for the first and second generation. 

3.5 Impact of the Sparticle Spectrum 

As already stated in the previous sections, the GUT scale values of the quark and lepton 
Yukawa couplings (and masses) strongly depend on the sparticle spectrum. We will now 
analyse its impact numerically in more detail. Because of the large number of relevant 
parameters, we do not attempt to discuss each of them separately, but rather make a 
parameter scan. 

For our scan, we take the three example ranges of SUSY parameters used for our ana- 
lytical estimates in section 3.1, which are listed explicitly in table 1. In addition, for tan/3 
we assume a range from 30 to 50. The sfermion mass parameters and At are scanned with 
a step size of 1 TeV (including At = 0) , the mass of the lightest gaugino was changed with 
a step size of 0.5 TeV and tan (3 is scanned with a step size of 10. Although this seems to be 
a rather coarse scan, we note that we have shown in section 3.1 that the extremal values of 
the SUSY threshold corrections correspond to the extremal values of the SUSY parameters. 
Therefore, increasing the number of plotted points would not lead to enlarged ranges for 
the GUT scale quantities. The parameter points for which y^ becomes non-perturbative 
have been dropped from our analysis. We note that for simplicity we have introduced a 
slightly more restrictive cut and included only parameter points where y^ < 1 at MsusY 
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(which ensures perturbativity up to Mqut but eventually removes a few allowed param- 
eter points with large but still perturbative yb)- Furthermore we have chosen MsusY to 
be 1 TeV. The masses of the first two sfermion generations have been assumed identical, 
which is inspired by universal high scale boundary conditions for sfermions. We note that 
a small mass splitting does not change the conclusions from this plot; however a large mass 
splitting can reduce the threshold effects due to a reduced value of the function | H2 \ in the 
formulae (2.2) to (2.8). 

The results of our parameter scans are presented as scatter plots in figure 5. The grey 
areas correspond to the (la) quark mass errors for each shown data point. In the first 
column, nis and ma have been varied, and in the second column yb and yt, using best- 
fit values for the remaining fermion masses. In comparison to the quark mass errors, the 
charged lepton mass errors are negligible. In figure 5 only those data points are shown 
where yi,, including its la error, stays perturbative. For comparison we have also included 
the best-fit values, which would be obtained without SUSY threshold effects. It can be 
seen that, with SUSY threshold corrections included, the shown ratios are (for almost all 
parameter points) significantly shifted to large values for case 5+ and to smaller values for 
case g^ and a. 

The plots also reveal some interesting features, which we discuss now. For example 
we notice that small values for m^/ms are correlated (nearly linearly) to small values of 
me/rrid. This correlation is connected to our assumption that the sfermion masses of the 
first two families are assumed to be identical. Because of the quark mass errors, this 
correlation is somewhat smeared. Looking next at the third family relations yr/Ub versus 
VtlVh^ we find that there is an additional tan/3 dependence, which is due to the fact that the 
relation of top mass and Yukawa coupling differs from the relations for down-type quarks 
and charged leptons by a factor of tan /3. For all three cases we can therefore distinguish 
three bands, which correspond to the three values of tan/3 in our parameter scan. The 
scans show that for case g- and a it is in principle possible to obtain third family Yukawa 
unification for tan/3 ~ 50 (in fact for tan/3 somewhat below 50 in case a), whereas for 
g^ we found that it could not be exactly realised. Although y^ ~ yt can be achieved 
for tan/3 = 50, At = —1 TeV, [i = 0.5 TeV, light gaugino masses, mj, = 1.5 TeV and 
wiq = rauo, = 0.5 TeV, we found yr/yb ^ 1-1 in the considered parameter range. 

We note that since we have not specified the remaining SUSY parameters (which do 
not enter the formulae for the threshold corrections) we have not applied various relevant 
phenomenological constraints on the spectrum. We would therefore like to warn the reader 
that for additional assumptions for the SUSY parameters some of the considered parameter 
points may be phenomenologically challenged already by the present experimental data (e.g. 
from B physics, g^ — 2 or from LFV charged lepton decays in seesaw scenarios). In particular 
the cases with large tan/3, small pseudoscalar Higgs mass and large \At\ may (under some 
conditions) be challenged by Bs -^ /U^/U~ data and the case /U < (and M2 > 0) may seem 
already disfavoured by ^f^ — 2 if the assumption is made that the (> 3a) deviation from the 
SM prediction is restored by SUSY loop effects. 

The ranges for the quark and lepton Yukawa couplings at the GUT scale and for the 
ratios of interest are presented in table 5 without quark mass errors. They can be compared 
to the results of our semianalytical treatment given in table 4. Comparing the two tables one 
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Figure 5: Scatter plots illustrating the ranges of GUT scale ratios jn^/rris, m^/md, yr/yb and yt/yb cor- 
responding to the example ranges of SUSY parameters in table 1 including At = Q (first row: case g+, 
second row: case g-, third row: case a) for tan/3 = 30, 40 and 50. The green dashed lines in the left plots 
correspond to the GJ relations. In the right plots they indicate the ratios j/r/yt. = 1 and yt/yb = 1. The 
green stars correspond to case (no SUSY threshold corrections). The point-bands in the plots on the right 
correspond to the different values of tan/3. The grey areas correspond to the (Icr) quark mass errors for each 
shown data point. 
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Case 


Case 5+ 


Case g- 


Case a 


me/md 


0.39 


[0.45 


0.55] 


[0.26, 0.34] 


[0.24, 0.36] 


m^/ms 


4.35 


[4.95 


6.03] 


[2.92, 3.80] 


[2.65, 3.92] 


Vr/yb 


1.32 


[1.23 


2.23] 


< 1.40 


< 1.35 


yt/vb 


1.93 


[1.67 


3.03] 


<2.17 


< 1.76 


ye in 10-^ 


0.88 


[0.87 


0.97] 


[0.80, 1.09] 


[0.92, 1.44] 


y^ in 10-2 


1.85 


[1.83 


2.06] 


[1.68, 2.31] 


[1.95, 3.05] 


Vr 


0.34 


[0.34 


0.39] 


[0.31, 0.44] 


[0.36, 0.60] 


yd in 10-4 


2.22 


[1.66 


2.10] 


[2.39, 3.80] 


[2.73, 5.73] 


ys in 10-2 


0.43 


[0.32 


0.40] 


[0.46, 0.73] 


[0.52, 1.10] 


yb 


0.26 


[0.16 


0.30] 


>0.23 


>0.29 


yu in 10-6 


2.75 


[2.73 


2.76] 


[2.74, 2.83] 


[2.75, 2.87] 


yc in 10-^ 


1.34 


[1.33 


1.35] 


[1.34, 1.38] 


[1.34, 1.40] 


yt 


0.50 


[0.48 


0.51] 


[0.50, 0.58] 


[0.50, 0.62] 



Table 5: Ranges for the GUT scale ratios and Yukawa couplings (corresponding to the example ranges of 
SUSY parameters in table 1, including additionally At — 0) from our numerical analysis with tan/3 — 40. 
Case refers to the case without SUSY threshold corrections. Quark mass errors are not yet included. The 
results can be compared with the semianalytic estimates of table 4. Our final results, which include the 
experimental mass errors, are given in tables 6 and 7. Where yt becomes nonperturbatively large we have 
given only the boundary, for which yt stays perturbative up to the GUT scale. 



can see that the "mean values" of the ranges agree well for the first two generations, however 
the extremal values are somewhat different. This is no surprise since in our naive estimates 
we have generically overestimated the ranges of the e^ (since we have ignored possible 
correlations between the corrections). For the third generation, the ranges (boundaries) are 
slightly shifted. This effect is caused by the modified RG running with SUSY threshold 
corrections included. For the first two generations this effect is smaller due to the smaller 
Yukawa couplings. 

While the quark mass errors are still ignored in tables 4 and 5, they are included in 
our final results listed in tables 6 and 7. Here we have varied all the quark mass errors 
simultaneously. Comparing case (without SUSY threshold corrections) to the cases g± 
and a, we see that the ranges for all types of Yukawa couplings, for down-type quarks, 
charged leptons and up- type quarks, as well as for all three generations, can be significantly 
affected by the SUSY threshold corrections. We note that for the up- type quarks, the 
changes are indirect in the sense that they are induced by modified RG running (mainly) 
due to the corrected fe-quark and r-lepton Yukawa couplings. 



4 Implications for GUT Model Building 

In the previous sections we have seen that comparatively wide ranges of Yukawa couplings 
are allowed at the GUT scale, if possible SUSY threshold corrections are taken into account. 
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From todays perspective, since we do not have any experimental confirmation for low-energy 
SUSY and therefore no knowledge of the sparticle spectrum, this can have a large impact 
for GUT model building. For instance, as has been shown in [15], the threshold corrections 
can make the GJ relations consistent with the latest experimental data on quark masses. If 
the GJ relations are assumed at high energies, this can be understood as a constraint on the 
SUSY parameter space and points to scenarios with M3 < and /x > (if phenomenological 
consistency with experimental results on (g — 2)^ is assumed to be restored by SUSY loops 
[25]). Another interesting aspect is that the wide allowed ranges for the Yukawa couplings 
at the GUT scale open up new possibilities for constructing SUSY GUT models to address 
the flavour problem. 

One example for an application of such alternative GUT scale ratios m^/iris and me/rnd 
can be found in [26] , where an approach has been presented to realise the phenomenologically 
successful relation 9i2 + ^c ~ vr/4 (so-called quark-lepton complementarity [27]) in unified 
theories. In this approach, the Yukawa matrices for the charged leptons and down-type 
quarks emerge from the identical higher-dimensional operators where quarks and leptons 
are unified in representations of the Pati-Salam gauge group G4221 = SU(4)cX SU(2)lX 
SU(2)rX U(1)b-l- After spontaneous breaking of 6*4221 to the SM gauge group Clebsch- 
Gordan factors lead to different (GUT scale) values for the charged lepton and down-type 
quark Yukawa couplings. In the approach of [26], for example, m^/nis = 2 was postulated 
at Mgut- 

More generally, the assumption that the Yukawa matrices for the charged leptons and 
down-type quarks are generated from the same set of higher-dimensional operators in quark- 
lepton unified theories leads to a large variety of possible ratios m^/ms and me/md which 
correspond to different choices of operators and their associated Clebsch factors. A table 
with a collection of possible Clebsch factors in the context of Pati-Salam theories can be 
found in the appendix of [28] . Any of these combinations of Clebsch factors which results 
in ratios m^/nis and ra^jrad consistent with the ranges in table 6 are a priori interesting 
new options for GUT model building. On the other hand, values of m^/ms larger than 3 
seem to be difficult to realise within the Pati-Salam framework, which may point to models 
where the leading contributions to charged lepton and down-type quark Yukawa matrices 
emerge from different operators. 

Finally, if low-energy SUSY is found at the LHC and if the SUSY parameters are 
determined at the LHC (and/or ILC), the SUSY threshold corrections can be calculated 
and significantly more precise statements about the GUT scale values of the quark and 
lepton masses can be achieved. In addition, more precise values of the low-energy quark 
masses (in particular of the first two generations) would be highly desirable to improve the 
GUT scale predictions. This would allow to select among the possible GUT scale ratios of 
quark and lepton Yukawa couplings and may point to unexpected new relations. 



5 Summary and Conclusions 

In this study we have investigated the effect of SUSY threshold corrections on the values of 
the running quark and charged lepton Yukawa couplings at the GUT scale within the MSSM 
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Ratio 


CaseO 


Case g^ 


Case g- 


Case a 


tan /3 = 30 


nie/md 


[0.28, 0.67] 


[0.30 


0.86] 


[0.21, 0.61] 


[0.20, 0.62] 


m^/nis 


[3.39, 6.07 


3.73 


7.79 


2.54, 5.49 


[2.40, 5.63 


Vr/yb 


[1.27, 1.38 


1.20 


2.02 


0.71, 1.43 


[0.60, 1.39 


yt/Vb 


[2.56, 3.02 


2.36 


4.19 


1.50, 3.28 


[1.14, 2.87 


tan (3 = 40 


nie/md 


[0.28, 0.67] 


[0.31 


0.93] 


[0.19, 0.59] 


[0.17, 0.60] 


mfj,/ms 


[3.39, 6.07] 


[3.85 


8.41] 


[2.28, 5.30] 


[2.07, 5.47] 


yr/yb 


[1.26, 1.38 


1.16 


2.32 


< 1.46 


< 1.41 


yt/yb 


[1.77, 2.11 


1.55 


3.31 


<2.38 


< 1.94 


tan /3 = 50 


rrie/md 


[0.28, 0.67] 


[0.32 


1.00] 


[0.16, 0.57] 


[0.14, 0.59] 


m^/nis 


[3.39, 6.07 


3.98 


9.06 


2.02, 5.12 


[1.72, 5.31 


Vr/yb 


[1.25, 1.38 


1.08 


2.73 


< 1.49 


< 1.43 


yt/vb 


[1.22, 1.50 


0.94 


2.74 


< 1.81 


< 1.31 



Table 6: Ranges for the GUT scale ratios nifj^/ms, rrie/md, yr/yb and yt/Vb corresponding to the example 
ranges of SUSY parameters g+, g- and a defined in table 1 including At — 0. The results have been extracted 
from the numerical analysis (parameter scan) with tan/3 = 30,40 and 50, where in addition to the SUSY 
threshold corrections the present experimental errors for the quark masses have been included. Case refers 
to the case without SUSY threshold corrections. Where yb becomes nonperturbatively large we have given 
only the boundary, for which yb stays perturbative up to the GUT scale. 



in the large tan f3 regime. We have therefore solved the set of RGEs from the top mass scale 
rrit to the GUT scale MguT) taking into account the one-loop threshold corrections at a 
SUSY scale (Mgusv) which has been assumed to be about 1 TeV. With superparticle masses 
not too widely split around MsusY) RG running can be performed in a two-step procedure, 
first with the SM RGEs from mt to MsusYj where one-loop matching is performed, and 
then from MgusY to Mqut using MSSM RGEs. The additional right-handed neutrino 
thresholds within the seesaw framework have also been considered but turned out to have 
no significant effects. 

At the matching scale MsusY) we have included all relevant tan /^-enhanced one-loop 
corrections to the quark and lepton Yukawa couplings, which stem from SUSY QCD contri- 
butions (for all down-type quarks), SUSY EW contributions (for all down-type quarks and 
charged leptons) and from the trilinear soft SUSY breaking coupUng At (for the b quark). 
For this purpose we have also calculated the SUSY EW corrections to the charged lepton 
Yukawa couplings in the EW unbroken phase. In general they cannot be neglected in a 
quantitative analysis of the SUSY threshold effects as we have shown. 

Using the above described approach, we have analysed which values of quark and charged 
lepton masses may be realised at Mgut- Since the possible values depend on a large number 
of SUSY parameters as well as on the quark mass errors, we have organised our analysis as 
follows: 

• We have first analysed the size of the SUSY threshold corrections for three example 
ranges of SUSY parameters g+, g- and a (listed in table 1). While there is a direct 
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Yukawa 


CaseO 


Case 5+ 


Case (7_ 


Case a 


tan/3 = 30 


Ve 


n 10"^ 


0.62 


[0.62 


0.67] 


[0.58, 0.66] 


[0.63, 0.79] 


y,. 


m 10"^ 


[1.30, 1.32 


[1.32 


1.41 


[1.22, 1.40 


1.34, 1.66 




yr 


0.23 


[0.23 


0.25 


[0.22, 0.25 


0.24, 0.30 


yd 


m 10""^ 


[0.92 


2.26] 


[0.75 


2.14] 


[0.98, 3.11] 


[1.06, 3.89] 


ys 


n 10"^ 


[0.21 


0.39 


[0.17 


0.37 


[0.23, 0.53 


0.25, 0.67 




yb 


[0.17 


0.18 


[0.12 


0.20 


[0.16, 0.34 


0.18, 0.48 


Vu 


m 10-^ 


[1.79 


3.88 


[1.78 


3.89 


[1.78, 3.93 


1.79, 3.98 


yc ] 


n 10"^ 


[1.14 


1.54] 


[1.13 


1.54] 


[1.14, 1.56] 


[1.14, 1.58] 




yt 


[0.46 


0.51 


[0.46 


0.52 


[0.46, 0.54 


0.46, 0.57 


tan/3 = 40 


ye 


n 10"^ 


[0.87 


0.88] 


[0.86 


0.99] 


[0.79, 1.62] 


[0.91, 1.77] 


y,. 


m 10-^ 


[1.83 


1.87 


[1.82 


2.08 


[1.67, 3.43 


1.93, 3.74 




2/r 


[0.34 


0.35 


[0.34 


0.39 


[0.30, 0.67 


0.36, 0.76 


yd 


m 10""^ 


[1.30 


3.21] 


[0.97 


3.05] 


[1.40, 8.41] 


[1.59, 9.81] 


ys 


n 10"^ 


[0.30 


0.55 


[0.23 


0.52 


[0.33, 1.44 


0.37, 1.69 




yb 


[0.25 


0.27 


[0.15 


0.33 


>0.22 


>0.27 


yu 


m 10-" 


[1.79 


3.91 


[1.78 


3.92 


[1.79, 4.07 


1.80, 4.15 


yc^ 


n 10"^ 


[1.14 


1.55] 


[1.14 


1.56] 


[1.14, 1.62] 


[1.14, 1.65] 




yt 


[0.47 


0.53 


[0.46 


0.54 


[0.47, 0.65 


0.48, 0.72 


tan/3 = 50 


ye 


n 10"^ 


[1.18 


1.23] 


[1.13 


1.53] 


[1.04, 2.31] 


[1.30, 3.80] 


yu 


m 10"^ 


[2.50 


2.60 


[2.39 


3.24 


[2.19, 4.89 


2.74, 8.04 




yr 


[0.50 


0.52 


[0.47 


0.69 


[0.42, 1.07 


0.56, 2.20 


yd 


m 10-" 


[1.77 


4.46] 


[1.20 


4.56] 


[1.91, 12.22] 


[2.36, 22.68] 


ys 


n 10"^ 


[0.41 


0.77 


[0.28 


0.78 


[0.44, 2.10 


0.55, 3.90 




yb 


[0.36 


0.42 


[0.19 


0.60 


>0.30 


>0.43 


yu 


m 10-" 


[1.81 


3.95 


[1.79 


4.00 


[1.80, 4.18 


1.81,4.21 


yc in 10"^ 


[1.15 


1.57] 


[1.14 


1.59] 


[1.15, 1.66] 


[1.16, 1.67] 


yt 


[0.49 


0.56 


[0.46 


0.59 


[0.48, 0.78 


0.50, 0.86 



Table 7: Ranges for the GUT scale values of the Yukawa couplings corresponding to the example ranges 
of SUSY parameters g+, g- and a defined in table 1 including At — 0. The results have been extracted 
from the numerical analysis (parameter scan) with tan/3 — 30,40 and 50, where in addition to the SUSY 
threshold corrections the present experimental errors for the quark masses have been included. Case refers 
to the case without SUSY threshold corrections. Where yt becomes nonperturbatively large we have given 
only the boundary, for which yt, stays perturbative up to the GUT scale. 
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dependence on the gaugino masses and on /i in the corresponding contributions in 
equations (2.2) to (2.8), one can see that the dependence on the squark and slepton 
masses is somewhat hidden in the function H2. We have used the ranges for the 
threshold corrections to estimate the possible variety of GUT scale values of the 
Yukawa couplings due to the unknown SUSY spectrum semianalytically. 

• We have next used best-fit values for the low-energy quark and lepton masses and 
isolated the effects of tan/3 and |/i| as well as of the signs of /U and of the gluino mass 
M3 (keeping Mi and M2 positive). On the one hand this illustrates the dependence 
on tan/3 and \^\ anticipated from the leading SUSY QCD contribution in equation 
(2.2). On the other hand, comparing the cases n > 0, M3 < to ^ < 0, M3 > (with 
sparticle masses fixed at this stage at 1 TeV) demonstrates the importance of the EW 
contributions to the threshold corrections. Our results are shown as contour plots in 
figures 3 and 4. 

• We have then discussed the effects of At, which can provide an important contri- 
bution to the correction for the b quark Yukawa coupling yh (but can be neglected 
for the others). Furthermore, varying MsusY (but keeping it still around the TeV 
scale) did in general not have a significant effect. Similarly, including right-handed 
neutrino thresholds, which appear in seesaw scenarios for small neutrinos masses, did 
not significantly affect the GUT scale values of the quark and charged lepton Yukawa 
couplings, at least compared to the relevant SUSY parameters. 

• Finally we have turned to the numerical discussion of the impact of the sparticle 
spectrum using the same example ranges of SUSY parameters as for the semianalytic 
estimates. We note that since we have not specified the remaining SUSY parameters 
(which do not enter the formulae for the threshold corrections) we have not applied 
various relevant phenomenological constraints on the spectrum (cf. discussion in sec- 
tion 3.5). The numerical results for our example ranges of SUSY parameters g^, 
g- and a are summarised as scatter plots in figure 5 for tan/3 = 30, 40 and 50 and 
additionally in table 5 for tan/3 = 40. 

Our final results, with quark errors included, are presented in tables 6 and 7. They 
provide the resulting possible ranges for quark and charged lepton Yukawa couplings at 
the GUT scale in the presence of threshold corrections (for the example SUSY parameter 
ranges g^, </_ and a) as well as the resulting ranges for several relevant quark and lepton 
mass ratios. The tables illustrate that with SUSY threshold effects (and quark mass errors) 
included, a wide range of GUT scale values of down-type quark and charged lepton Yukawa 
couplings could be realised, consistent with the low-energy experimental data on quark and 
charged lepton masses. One interesting aspect is that this opens up new possibilities for 
constructing GUT models of fermion masses and mixings. SUSY threshold corrections can 
on the one hand improve consistency with existing postulated GUT scale relations such as 
the GJ relations m^/ms = 3, me/rrid = 1/3 or third family Yukawa unification yt = Vb = Vt, 
but they might also point to new GUT scale relations (cf. discussion in section 4), if low- 
energy supersymmetry is discovered at the LHC. 
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